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PREFACE 


The wo^k described herein was perforaed by the Control and Energy 
Conversion Division of the Jet Propulsion Laboratory. 
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DEFINITION OF SYlffiOLS 

concentration of gas A, g*mol/CB3 
drag coefficient 

concentration of A in the gas bubble, in the cloud-wake 
region, in the emulsion phase, and at the particle surface, 
respectively, g mol/cm^ 

average concentration of A in the bubble, g moI/cm3 

average concentration of A in the bed 

concentration of A in the entering gas stream, g mol/cm3 

molecular diffusion coefficient of gas, cm^/sec 

effective bubble diameter, cm 

particle diameter, cm 

surface mean particle size, cm 

tube diameter, cm 

feed rate of solids, outflow rate of solids, and carryover 
rate of solids by entrainment, respectively, gm/sec 

2 

Froude number at minimum fluidization condition, Ugjf/dpg, 
dimensionless 

980 cm/sec2, acceleration of gravity 
height, cm 

overall coefficient of gas Interchange between bubble and 
emulsion based on volume of bubbles, volume of emulsion, 
and total volume of bed, respectively, sec“^ 

entrainment coefficient 

mass transfer coefficient between fluid at a particle, cm/sec 

mass transfer coefficient in fluidized bed, sec~1 

overall mass transfer coefficient for bubbling bed, sec"^ 

mass transfer coefficient between fluid and a particle at 
minimum fluidization, cm/sec 

mass transfer coefficient for particles in bubbles, cm sec”^ 

iv 
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i height, OB 


Lf 

Lnf 

Mw 

1 

H 

NAb 

AP 

Pq. ’1. 
P2. Pb 

*'o 

fp, Bihf 

Tp, maxf 
R 


Rear 

Sc 

Sh 

Sht 


height of a bubbling i' .uidized bed, cm 

bed height at miniaum fluidization condition, cm 

molecular weight 

bubble frequency, sec~^ 

number of moles of A 

number of moles of A in bubble 

pressure drop, gm*wb/sm^ 

size distribution of feed solids, outflow solids, entrained 
solids, and solids in the bed, respectively, cm~^ 

particle radius, cm 

minimum particle radius in feed particles, cm 

maximum particle radius in feed particles, cm 

half the distance between particles in particulate bed 
or emulsion, cm 

particle Reynolds number, dpUQ^/M, dimensionless 

particle Reynolds number at minimum fluidization condition, 
dpUaf^g/^, dimensionless 

Schmidt-number , /t//»gO, dimensionless 

Sherwood number, k^jdp/D, dimensionless 

Sherwood number for particles in bubbles, k^^^dp/D, dimen- 
sionless 


Shgif Sherwood number for particles in emulsion or at minimum 

fluidization conditions, K^,nfdp/D, dimensionless 

Up velocity of a bubble rising through a bed, cm/ sec 

Upp velocity of a bubble with respect to the emulsion phase, cm/ sec 

Uf Inward velocity of gas at minimum fluidizing conditions, cm/sec 

Ugif superficial fluid velocity at minimum fluidizing conditions, 

cm/sec 


Up terminal velocity of a falling particle, cm/sec 


V 



77-25 


«0 


% 

W 

a 

a' 

rc 


superficial fluid velocity (measured on an empty tube basis) 
through a bed of solids, cm/ sec 

volume of a gas bubble, cm^ 

weight of solids, g 

specific surface or surface of solid per volume of 

specific surface or surface of solid per volume of 

ratio of solids dispersed in bubbles to the volume 
In the bed, dimensionless 

ratio of solids in the cloud-wake region to volume 
in the bed, dimensionless 

ratio of solids in the emulsion to volume of bubbles in 
the bed, dimensionless 


bed, cm“^ 
solid, cm*' 
of bubbles 

of bubbles 


8 fraction of fluidized bed consisting of bubbles, dimensionless 

€ void fraction, dimensionless 


^e* ^f> fraction in the emulsion phase of a bubbling bed, 

^mf> ^m ^ bubbling bed as a whole, in a bed at minimum fluidizing 

conditions, and in a packed bed, respectively, dimensionless 


Pgt ^ 
^s 


viscosity of gas, g/cm*sec. 

density of gas and solid, respectively, g/cm^ 

sphericity of a particle, dimensionless 


Vi 
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ABSTRACT 


Modeling of the fluidized bed for silicon deposition is described. 
The model is intended for use as a means of improving fluidized bed 
reactor design and for the formulation of the research program in 
support of the contracts of the Silicon Material Task for the development 
of the fluidized bed silicon deposition process. A computer program 
derived from the simple modeling is also described. Results of some 
sample calculations using the computer program are shown. 
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SECTION I 
INTRODUCTION 


The Low-Cost Silicon Solar Array (LSSA) project sponsored at 
JPL by ERDA has the responsibility for developing the technology and 
the industry -capability for producing low cost solar arrays. The objective 
in the achievement of a production capacity of 500 megawatts (peak) 
at a cost of $500 per kW by 1986. The program of LSSA is divided into 
four technology tasks, each of which has the responsibility for advancing 
a particular technical area so as to achieve the production and cost 
goals. The first of these tasks is the Silicon Material Task, for 
which the 1986 objective is to reduce the price of Si material, suitable 
for the fabrication of solar cells having satisfactory performance, 
from the present price of about $65 /kg (for semiconductor grade Si) 
to less than $10 /kg. 

To achieve the price objective of the Silicon Material Task, 
new areas of chemical processing and solar cell technology must be 
studied and developed, since no modification of the basic processes 
presently used commercially will lead to the attainment of the $10/kg 
price. The inability of the present commercial process to meet the 
task objective is revealed from a consideration of the features of 
this process. The first step is the reaction of metallurgical grade 
Si, having a purity of about 98.5}, with KCl in a fluidized bed reactor 
to yield the intermediate SiHCl 3 ® mixture of chlorinated silanes. 

After chemical treatment and fractional distillation to obtain very 
pure SiHClj, semiconductor grade Si .'.s produced by chemical vapor depo- 
sition resulting from the H 2 reduction of SiHCl 3 . This method of producing 
extremely pure Si is costly and energy intensive, the final deposition 
consuming an estimated 385 kWh/kg Si product. 

Several different chemical processes are being investigated under 
JPL contracts in the task program. Two of these, at Union Carbide 
and Battelle, incorporate fluidized bed reactor technology. The contract 
with Union Carbide involves the hydrogenation of SlClji to S 1 HC 13 in 
a fluidized bed as well as a process for the deposition of Si from 
SiHi| in another fluidized bed. In the contract at Battelle the process 
is for the Zn reduction of SiCl 4 in a fluidized bed reactor. These 
two efforts are presently the most advanced in the task program and 
are being appropriately emphasized. 

To support the developments in fluidized bed technology, which 
is the basic element of these two important contracts, the task program 
includes an in-house JPL subtask for studies of fluidized bed reactor 
technology as well as a subtask for consultation in this area by Professors 
0. Levenspiel and T. Fitzgerald of Oregon State University. The in- 
house program consists of efforts in modeling of the deposition of Si 
in a fluidized bed reactor, experimental studies of SiHi| pyrolysis 
and fine particle fluidization, and thermodynamic and chemical engineering 
analysis of the reactions in the fluidized process for silicon production. 
The initial phase of the effort in modeling is described in this report. 
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SECTION II 
APPROACH 


The purpose of this study is to develop an analytical model of a 
particular chemical reaction — the pyrolysis of SiHjj— occurring in a 
fluidized bed reactor. The model is intended for use in determining 
guidelines for reactor design and for the structuring of experimental 
investigations to achieve system optimization. The methodology of this 
study is (1) to specify the system, (2) to describe the particular and 
general assumptions, (3) to state the appropriate physical laws to be 
used at each step, (4) to apply the limitations of the system characteristics 
and assumptions in the derivations of the mathematical expressions, and 
(5) to develop a suitable computer program. To facilitate the systematic 
development of the computer program, this report is divided into sections 
for (1) a general description of the process, (2) modeling of the gross 
bed behavior, (3) modeling of the particulate bed, (4) modeling of the 
bubbling bed, (5) a description of gas interchange between bubble, cloud, 
and emulsion, (6) an analysis of mass transfer rate to solid particles in 
a bubbling bed, (7) particle growth rate in a particular bed, (8) particle 
growth rate in a bubbling bed, (9) entrainment, (10) overall mass balance 
of the bed, and (11) results, conclusions, and recommendations. 

This paper deals with the fluidized bed silicon deposition process 
designed to eliminate the difficulties associated with the current 
state-of-the-art. In Figure 2-1, a schematic diagram of a fluidized 
bed is shown. From the bottom left-hand side, small feed-particles 
of silicon are Introduced while the feed gas of silane and hydrogen 
(carrier gas) comes in from the bottom of the bed. The silane is pyrolyzed 
at the surface of the silicon particles, the resulting Si deposits 
causing the Increase in size of the seed particles, which are eventually 
removed as the product. Unreacted silane and hydrogen gas escape through 
the top of the bed concurrently elutriating small amounts of srnal^- 
size particles. This process is continuous and has the advantage of 
using a very large reaction surface area to produce pure silicon on 
a much larger scale than the conventional chemical vapor deposition 
(Siemens) process. 

Our problem can be stated as follows: Given a set of operating 

conditions, which Include the bed size, the feed particle characteristics, 
and the feed gas characteristics, determine (1) the type of fluidization 
(l.e., whether it is bubbling or particulate type), (2) the optimum 
gas flow rate, (3) the power requirement, (4) the silicon particle 
growth rate, (5) the entrainment rate, and (6) the deposition rate 
of Si. 

Several assumptions and a number of limiting conditions were employed 
in the modeling of this process. They are as follows: 
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Figure 2-1. Fluidized Bed Concept for Tilicon Production 


(1) At temperatures above 1000® C, the primary reaction is the heter- 
ogenous pyrolysis of SiHn on the surface of the silicon particles. 

The rates of adsorption and reaction are much faster than 
the transport of the reactant species, SiH^, to the reaction 
zone (the hot surface), and thus the growth rate of the 
particle size is diffusion-limited. 

(2; The bed is either bubbling or partlculately fluidized. In 
other words, the fluidization of particles smaller than 50 
which would cause abnormal fluidization, will not be considered. 

These two types of fluidization are schematically shown in 
Figure 2-2. For larger particles, bubbles moving up are 
observed (Figure 2-2a), and for i^maller particles (but still 
larger than 50 /*) uniform fluidization is observed (Figure 
2-2b). This uniform phase in the particulate bed and the 
bubbling bed (outside the bubbles) is called the emulsion. 

The basic flow chart for the modeling is given in Figure 2-3. As shown 
here, first, the set of input parameters is given. From this, the gross be- 
havior of the bed is obtained including the operating bed weight, the pressure 
drop across the bed, the distributor requirement and the type of fluidization. 
If the bed is particulate fluidized, the mass transfer rate of gas onto solid 
particles and thus the particle growth rate can then be obtained relatively 
simply. If the bed is bubbling, the bubble and the emulsion characteristics 
must be studied separately with the aid of pertinent experimental data. Then 
the interaction between the two phases must be Investigated to determine the 
mass transfer rate of the gas onto the solid particles in order to permit the 
calculation of the particle growth rate as a function of the particle diameter. 
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tUiniNG PaR rulATE 



Figure 2-2. Types of Fi idization 


When this growth rate Is obtained, the mass balance of silicon 
particles in the feed, product, entrainment streams and the bed is 
then calculated. The partible size distribution of the bed and the 
deposition rate can also be determined. 



Figure 2-3. Flowchart for Modeling 
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SECTION III 

MODELING OF GROSS BED BEHAVIOR 


The input parametera necessary for the definition of the fluidized 
bed are as follows: feed particle size distribution pQ(dp); production 

rate F^; sphericity of particles which is assuaed to be unifora 
regardless oi the particle size; void fraction in packed bed cb, as 
a function of ^3 (Ref. 3-1); voidage at ainiaua fluidizing condition 
cgf, and as a function of ^3; gas density p.; solid densitv pg; viscosity 
p; diffusion coefficient 0; concentration of silane in the feed gas 
Cgi; bed weight W; and the bed diameter dt> 

In order to estiaate the bed characteristics, it is necessary 
to know the average particle size in bed. However, without the Know- 
ledge of the bed particle size distribution P^Cdp) this is not possible. 
Thus, at the beginning, 7p is obtained froa Pg. This will be replaced 
by ^p(Pb) later. Since the important paraaeter in the bed is the ratio 
of surface area to volume of solid particles, which is inversely pro- 
portional to dp, the average particle size is obtained as the inverse 
of the average of the inverse of the particle size. That is. 


1 



Using the above inforaation, the gross behavior of the fluidized 
bed is then predicted as follows: 

(1) Bed height at mi.iifflua fluidizing conditions: There is 

a certain velocity of gas for the given bed below which 
the solid particles are pac'ted together and no fluidization 
is observed. Above this velocity, the particles are moving 
around in the bed with gas to create a fluid-like phase 
in the bed. This velocity is called the minimuai fluidizing 
velocity, and the requirerent to reach this velocity is 
called the minimum fluidizing condition. 


W 4 1 


( 2 ) 


where g is the gravitational constant. 

(2) ihe pressure drop across the bed is obtained by balancing 
the pressure loss with the bed weight. 
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Ap s 



(3) 


(3) TIm superficial fluid velocity at ainiaua fluidizing condi- 
tions Ugf can be obtained by balancing the pressure loss 
Mitb either the viscous force for the saall-partiele»size 
bed or the kinetic force for the large-part iele<>size bed. 
For an assuaed value for the superficial fluid velocity 
of the bed Uq, then 


Oaf * 


V (/>s - /»g)g 

1650^ 


when Rep 


^p Pg, ^0 

Ji 


< 20 


(4) 


or 


2 5p (.Pa - Pg)g 

0inf = — idjen Ren > 1000 (5) 

a.5 


If 20 < Re < 1000, an Interpolated value of would be 
used. 


The superficial fluid velocity of the bed Uq can be esti- 
aated from this information. If <> vigorously mixing bed 
is assumed for maximum silicon production, then, according 
to literature observations (Reference 3*1): 


Uo > 2«mf 


( 6 ) 


however, there is an upper limit fo.’ Uq, since the Uq is 
larger than the terminal velocity of the particles, then 


(Pa - ^g)^ 
^ ' \ 3?iCd / 


where dp is the minimum particle size in bed and C^] can 
be obtained empirically; then some particles are lost through 
entrainment. This tendency of entrainment would be increased 

as Uq increases. 
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(4) The type of fldldization of the bed can be determined by 
an empirical formula (References 3-2, 3-3) under normal 
fluidization conditions. If the product of four dimensionless 
groups 


^s " ^tf 


is larger than 100, the bed is bubbling; otherwise, the 
bed is particulately fluidized. Here 


Fr«f 



( 8 ) 


u^f . dp • />g 

Rep.f = 


(9) 
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SECTION IV 
PARTICULATE BED 


If the bed is particulately fluidized, which is depicted in Figure 
4»1 , the mass transfer rate of gas onto the solid particles is obtained 
in the following way (Ref. 4-1). The individual particles of radius 
rQ are considered to be separated froa other particles by a distance 
of 2R, where R can be derived easily from the voidage of the bed: 


R 


*•0 

(1 -€)1/3 


( 10 ) 


Gas enters froa the bottoa with the velocity Uq. The aaount of gas 
transferred to the individual particle surfaces by diffusion, as well 
as by this forced convection, is to be calculated. 

This is a very coaplex two-dimensional flow problea for each 
particle. The description can be approximated by use of the penetration 
theory as follows. Assume that at t s 0, the concentration of the 
silane at r = R is the sane as that of the incoming gas, that at 
r : rQ, the concentration is zero, and that the diffusion irtiich is 
one-diaensional in the radial direction only proceeds until t s T. 

At t s T, the original situation is reestablished and the ensuing 
diffusion process is repeated aigain. This is repeated at everv bT 
interval where a is an Integer. The period T is proportions to the 
ratio between the incoalng velocity Uq divided by particle rauius rQ. 

The result is sunnarlzed in a formula for the mass transfer 
coefficient, k^j, 


^ N 

f \ 



Figure 4-1. Particulate Bed 
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Sh s 


“d” 


r 242 (1 - *1 

2 ♦ I ... — ■ — ■ ■ ■ 2 I tanh ^ 

Lt1 - (1 - «f)^/3]2 J 

— - tanh 4 

1 - (1 - «f)^^3 


( 11 ) 


where 


and 


1 1/2 1/2 

4 s rr : - ^ X 0.3 Rep Sc (12) 

(1 - «f)^/3 


))OPg 

"*p * — s — 


Sc 


PD 


At the limit when R*— 00, 


1/2 1/3 , , 

Sh^2 ♦ 0.6 Rep Sc (13) 

which is known as a single-sphere diffusion formula. This result will 
be used liter in the bubble phase analysis. 
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SECTION V 
LUBBLINQ BED 


To obtain the mass transfer rate of gas in the bubbling phase 
it is necessa.y to understand the bubble and the emulsion characteristics 
first. From the literature (Reference 3-D it is noted that the emulsion 
phase behaves very much like a bed at minimum fluidizing condition. 

Thus, the upward velocity of gas in the oaulsion phase is defined as 
Uf. Then, by assuming that the emulsion phase behaves approximately 
according to the minimum fluidizing condition, one obtains. 


Uf 


Umf 

*mf 


(1M) 


The. bubble phase is analyzed as follows. First the bubble diameter 
d^ is defined assuming that d{, is constant throughout the bed (although 
this is not true in the beginning of the bed), and that the bubble 
frequency passing by a certain point in the bed is n. Then the velocity 
of the rise of bubbles Ut)|. is similar to the velocity observed in the 
system of a bubbling liquid of low viscosity, in idiich case inertial 
forces predominate. By equating the inertial force to the gravity 
force, the relation for a single bubble is obtained: 


Obp = 0.711 (gdb)^/2 


(15) 


For bubbles rising together in the bed, the absolute velocity would 
be, then, 

Ub * <«0 - «mf> + Ubr- 

Depending on the bubble velocity there are two types of bubble 
formations as shown in (Figure 5-1 )• 

If Ub > Uf, a cloud will form around the bubble (Figure 5-1a) 
in which the gas circulates and never escapes to the outside. The 
thickness of the cloud decreases as Ub becomes larger. For vigorously 
bubbling beds, where Ub > 5Uf, the cloud thickness becomes negligible. 

If Ub < Uf (Figure 5- 1b), no cloud forms and the gas passes through 
the bubble as a shortcut to the top of the bed. The higher mass transfer 
rate is obtainable for the first case; i.e., Ub > Uf and this condition 
will be used. 

The superficial velocity Uq can be expressed as follows: 


Uo = (1 - 8) Unf + SUb 


(17) 
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CLOUO 


lUWLE 



Figure 5-1 . Bubbles 


trtiere S the voluae fraction of bubbles in the bed. The first tern 
is the velocity contribution from the emulsion phase and the second 
is from the bubble phase. 

Rewriting the above equation (17)» one obtains, for vigorously 
bubbling beds: 


% = 


Uq - (1 - 8)Ug,f _ Uq - Umf 

8 ^ 8 


( 18 ) 


If the height between two successive bubbles is h, then 


n = 


% 

h 


(19) 


and 


- 3 

^ 6<lb 2 db 



( 20 ) 
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The values of d^, (I5, 8, and h can be determined from (15), 

(16), (17), (19), and (20), with a known value of n. The value of 
n can be obtained from experimental data. 

On the other hand, the tendency for the bubbles to break up In- 
creases with bubble size. This limiting size Is computed as follows: 
idien the bubble breaks up. It does so In such a way that the solid 
particles go up through the center of the bubble and cause the separation 
of one cavity Into two smaller cavities, as shown In Figure 5-2. 

The upward velocity of gas throu^ a large rising bubble Is assumed 
to be nearly U^p. If U^p > U|^, then small entrained particles go up 
throu^ the middle of the bubble and the bubble breaks up. For this 
condition, a new n should be Introduced to make Ut,p = 


Figure 5-2. Mechanism of Bubble Breakup 
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SECTION VI 

GAS INTERCHANGE BETWEEN BUBBLE, CLOUD, AND EMULSION 


The flow pattern of gas through the fluidized bed can now be 
obtained from the above analysis. First, mass transfer coefficients 

^^be^bi (^bc)bf (^ce)b defined in the following way, which 
includes consideration of diffusion and gas circulation: 


1 dNAb ^ 

z ~r ~ '•^be)b • (CAb - CAe) 

''b dt 

= ^*^bc)b • (CAb - Cac) 


~ ^*^ce)b • (Cac * CAe) 


( 21 ) 


Here is the volume of bubble; N^^v is the number of moles of silane 
in the bubble; and CAb, Cac, and CAe are the mean c^ centrations of 
silane in bubble, cloud, and emulsion, respectively. Thus (K|^)(, is 
the volume of gas going from bubble to emulsion divided by volume of 
bubbles in bed per unit time. The other two are defined similarly. 
From the above equations: 


1 1 1 
-i — ■ — S ■ i ' + — — — 

^^be)b (Rbc)b (Kce)b 


( 22 ) 


(Rbc)b now obtained in the following way; There are two modes 
of gas movement between the bubble and the cloud. The first is gas 
flow into and out of a single bubble, and the second is the mass transfer 
by diffusion. A&*.uming a Davidson's model (Reference 3-1) for bubbles, 
then 


"mf 


(*bc)b s 4.5 5.85 


d 1/2 g1/4 

d5/4 
b 


and 


(*^ce)b * 6.78 


1/2 
<mfD0b 

db 


where (Kq*)j, involves the diffusion process only. 


(23) 


( 24 ) 
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SECTION VII 

MASS TRANSFER RATE TO SOLID PARTICLES IN BUBBLING BED 


The mass transfer rate to solid particles can now be derived 
using yt) as the volume of solids dispersed in bubbles divided by volume 
of bubbles. The terms emd y^ are similarly defined in the Definition 

of Symbols. The quantity yt is experimentally determined. 

For convenience the mass transfer process is considered to be 
an absorption by the solid of the silane gas present in the fluidizing 
gas stream. Since the flow of gas in the emulsion is very small, its 
minor contribution to total flow can reasonably be ignored. Thus, 
the fresh gas enters the bed only in the form of gas bubbles; moreover, 
for steady-state operations the measure of adsorption of silane is 
given by the decrease in its concentration within the bubbles. The mass 
transfer coefficient for such processes can be reported in a number 
of ways. 

Based on total surface of particles, an overall mass transfer 
coefficient (cm/sec) can be defined as 


1 doAb 
surface dt 


~ ^*^d^b ^^Ab ” ^As^ 


(25) 


or in terms of the falling concentration of silane in the rising bubbles 


dCAb dCAb surface of solids 

dt di ^*^d^b volume of bubble 


( 26 ) 


— g— <C*b - C*3) = (Kd)^, (CAb - Cas) 


where CAb mean concentration of silane in the bubble, Cas is 

the proper concentration measure of silane at the surface of the solid, 
and a s 6(1 - Cf)/d/^g. 

' The morphology of the silane present in a bubble as the bubble 
rises through the bed is considered now. Some of the silane is adsorbed 
by solids within the bubble; another portion is transferred to the 


^Although expressed in bubble-related terms, this equation covers the 
deposition process in all phases. 
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cloud where a part of it is adsorbed and where the rest is transferred 
further into the emulsion. The silane in the emulsion is assumed to 
be completely adsorbed, which is reasonable in light of the long contact 
time in the emulsion and the rapidity of the overall adsorption process. 
Thus we have 


adsorption in the transfer to 
overall adsorption = doud by solids ♦ the cloud 


where 

transfer to adsorption in the transfer to 
the cloud “ cloud by solids the emulsion 

and where 


transfer to adsorption in the 
the emulsion = emulsion by solids 

In symbols, the above becomes 


^^kb 

- — = -«b -^ = (Kd)b (CAb - Cas) (27) 


” 7b^d,t®^ (('Ab “ («As^ 
♦ ((*bc)b((-Ab “ (*Ac^ 

and 


(*(bc)b (('Ab ~ ('Ac) = yc*^d,mf®^ (Cac “ ^As^ * (*(bc^b (('Ac " ('Ae^ (2®) 


(^ce)b (('Ac ” ^Ae) “ ye*(d ,mf ( ('Ae “ ^ks^ 


(29) 


and 

surface area 6 

a' = = (30) 

volume of solid <|>gdp 

By eliminating C/^q and Cac from Equation (27), the general expression 
for the mass transfer coefficient, or Sherwood number, is obtained 
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(Kd> 


b 


- €f)D 

“ 'T ""^**overall 

8^s<*p 


* ^bBa 


Sht 

Shaf 


1 

. ^ 

(Kbc)^ 1 

b y^jBd + 

1 1 


+ 7e®<j 
b 


( 31 ) 


where 


Bd 


60 

P 


Shnf(3ec"^) 


(32) 


Thus with known values for Sh^ and Shg,f, and calculated gas Interchange 
rates, mass transfer coefficients and Sherwood numbers in fluidized 
beds can be estimated. 

K(] ^ is the mass transfer coefficient obtainable by using the 
single sphere diffusion formula 


*^d,t<*p 1/3 1/2 

Sht = — ^ — = 2.0 ♦ 0.6 S(j B«p,t 


(33) 


where 



B®p,t 


dpUt^g 

“*7 
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SECTION nil 

GROOTH RATE IN A PARTICULATE BED 


Based on total surface of particles, an overall mass transfer 
equation can be written 


1 dN^ . volume of bed 

surface of particles surface of particles dt 


^sdp dC/L 

dO - «f) dt 


= Kd(CA - Ca) 


Since Cg = o, due to the fact that this case is diffusion-controlled 
(Reference 8-1), the above equation can be easily integrated. The 
above equation can be rewritten 



Kd"CA 


where 


s K(j • g(1 - ^f)/^sdp 


dlnC|^ - -dl 


^A * ^Ai exp (IKd'l) 


'The average concentration of silane in the bed which could be used 
in average bed deposition calculation is 


C* s 


Cai( 1 - e'^'d Lf) 
^d 


(35) 


Therefore, the average mass transfer in the bed per unit surface 
of solids is given by 
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1 dNA 

- ' ■ ■ — - — * IcjiCa 
surface dt 


( 36 ) 


Finally! the particle growth rate Is calculated from this result as 


dR 

(It 


1 hw 

surface dt Pg 


^d^A 


(37) 


where Mw Is the molecular weight of silicon, and Pg Is its density. 
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SECTION IX 

PUITICLE QROVTH RATE IN A BU3BLING BED 

A siiiii&r procedure can be folloued to obtain the particle growth 
rate for the bubbling bed, except that ^*^d^K should be used in 
place of Do, k(j. The result is that the average concentration of silane 
in the bubbles is 


CAb 


% • 




( 38 ) 


The average aass transfer in the bed and the particle growth rate are 
given as, 


1 8N*b 
surface dt 


= O^d^K C*b, 


( 39 ) 


dR . Mw 

- • <kd>, . - 


( 40 ) 
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SECIION XI 

OVERALL MASS BALANCE OF THE BED 


The following definitions are used: 

Fq feed rate, F^: overflow rate, 

F2 elutriation rate 

Pq partic'e size distribution in feed 

Pi particle size distribution in overflow 

P2 particle size distribution in elutriation 

P5 particle size distribution in bed 


If Fq and Pq are known, then there are five unknowns to be cal- . 
culated and these are obtained by the following five relationships 
(refer to Reference 10-1): 

(1) Pj = P5, for a well-mixed fluidized bed 

(2) Ff F2 - Fq s total solid generation in bed. 

(3) From the definition of k, it follows that (see Reference 
10 - 1 ) 


P2(R) = k(R) HPb (R)/F2. 

(4) By definition, 



(5) In unit time, the mass balance on particles of size between 
R and R dR is given as follows: 


( solids 
entering 
feed 


\ / solids \ 

in j - I leaving in ] - 
J \ overflow / 


( solids growing out 
of the interval to 
a larger size 


( solids \ /solids growing into 

leaving in ] I the Interval from 
carryover J Va smaller size 


/ solids generation \ 
/ due to growth j 
\ within interval J 


\ 
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This rewrites as: 


- FlPl(R) - Wc(B)Pi(R) 


(R) P,(R) 3M dR 

- H ♦ — PWR) — 

dR R dt 


(R) s 0 


(41) 
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SECTION XII 
CALCULATION SAMPLES 


A computer program, incorporating the mathematical model described 
so far, was written. A listing of the program is attached as an appendix. 
The properties of the gah stream and silicon particles, the composition 
of the gas stream, bed diameter, initial bed weight, particle size 
distribution of feed particles, production rate, maximum particle size 
in the bed, the differential particle size interval and the number 
of orifices in the gas distributor were the input parameters to the 
model. The initial value for the t 'erage particle size in the bed 
is assumed to be the harmonic average of the feed particle size dl.«*--**'u- 
Lion. The computer program was found to converge to final value of 
average bed particle size in two iterations. The mat ^statical modeling 
equations were organized in various subroutines. The main program 
defines the variables, arrays and functions, sets the starting values 
for computation, provides the calling sequence for the subroutines, 
and outputs the calculated data. The output of the program consists 
of gross fluidization characteristics, bubble characteristics, mass 
transfer coefficients, particle size distributions and mass flow rates 
of feed, outflow or entrainment and overall growth rate of particles. 

The minimum particle size of the feed p^nin.f and the superficial velocity 
Uq are varied in three calculations as shown in the table for simplicity. 
The particle size distribution of the feed was assumed to be uniforik. 

Table 12-1 shows the results of some sample modeling calculations using 
this computer program. 

The results show that for the given set of values of variables, 
the minimum fluidizing velocity U„f, the terminal velocity U^, the 
bed height Lf, the necessary feed rate Fg, the entrainment rate F2, 
and the average particle size dp are listed. The convergence of the 
computation is very fast; at the second iteration, the asymototic 
solution is obtained. 

Fi was chosen after several trial and error calculations so that 
smooth bed operation is achieved. Such choice }f numbers is delicate 
due to the simplistic approach taken in the modeling. 

As compared to the first example, the second example shows the 
results when the superficial velocity is increased. This shows a smaller 
amount of Fq necessary than in the first case due to more deposition 
from the higher flow rate. The third example is for the higher average 
feed particle size, showing the reduced reaction due to smaller surface 
area. 


The results of the computer program were listed against experimental 
fluidization data obtained from fluidization of silicon particles with 
argon in a 1«ln. fluidization column. The results of the comparison 
are shown in Table 12-2. The mode of fluidization indicated by the 
model is consistently wrong. It shows that the simple Froude number 
criterion used in the model is not sufficient to predict t he mode of 
fluidization. The agreement between bed pressure drop data and calculated 
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Table 12-1. Examples of Caloulational Results 


e Input Paraaeters 


d^ = 5 cm 

^Ai 

= 5.0 X 10-® g-mole/cm3 

Pg = 1.23 X 10-** g/cm3 

*’p,nax,f 

= 75 m 

Pg = 2.4M g/cm3 

H 

= 200 g 

P = 3.*»5 X 10-** poise 

Fl 

= 0.92 g/sec 

D = 5.809 cm^/sec 




e Results 



‘*p,min,f = 25 m 

Oq = 19.2 cm/ sec 24.0 

25 M 
cm/ sec 

45 M 

19.2 cm/ sec 

Umf 

3.78 cffl/sec 

4.69 

5.36 

Ot 

9.62 cm/ sec 

9.62 

31.2 

Lf 

8.5 cm 

8.5 

9.4 

FO 

0.032 g/sec 

0.023 

0.048 

P2 

0.011 g/sec 

0.010 

0.009 


150.1 M 

167.2 

178.8 


bed pressure drop is excellent. The agreement between experimental 
and calculated values for minimum fluidization velocity is rough. 

This discrepancy may be due to several reasons. The model equation 
used for minimum fluidization velocity is a simplified low Reynolds 
number approximation of the more rigorous form. The experimental data 
itself was prone to considerable error. In many runs (1, 2, 3, and 7) 
the minimum fluidization velocity was determined by visual observation 
of bed behaviour. In cases where the UQ,f was determined from combined 
packed bed and fluidized bed data (6, 8, 9 and 10) the agreement is 
much better. Considering the possible significant experimental error 
in determining the minimum fluidization velocity, the computer model 
predicted the fluidized bed performance quite well. 
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SECTION XIII 
CONCLUSIONS 


As Dtentloned before, this modeling was intended for a simple 
description of the fluidized bed. Thus there are many assumptions 
and treatments which may not be practical. Further Improvements for 
the model will be conducted as experimental data become available. 

One such improvement area is the modeling of the fine-particle 
fluidized bed modeling. Although it may not be necessary if one wants 
to use larger-size particles when recycling some of the products by 
crushing them and using as the bed, it would still be beneficial to 
have this option of using fine particles. This modeling would be difficult 
at this time because very few experimental data are available. However, 
the JPL in-house experimental work is likely to yield the needed information 
for this important area of fluidized bed silicon production. 

Another improvement may come from consideration of gas interchange 
rate at the surface of solid particles. Unlike the implicit assumption 
that the same number of moles of gas come out as the number of moles 
of gas going in, two moles of hydrogen come out for each mole of silane. 

This might help the fluidization of fine particles, too. 

This report was reviewed by Professor T. Fitzgerald of OSU, acting 
as a consultant to LSSA. Dr. Fitzgerald (Reference 13-1) suggested 
that it might be useful to control the bubble size by artificial means: 
for example, vertical rods in the bed to break the bubbles larger than 
the size of the spacing. The reason is that, in the fluidized bed 
modeling, the prediction of the size of the bubble is the most difficult 
task (the error could be as high as 500%). This would improve the 
accuracy of the modeling. 

He also suggested that the following assumptions may prove trouble- 
some, and may be modified to fit real data, althou^ these considerations 
will make the model complex, and the improved accuracy may or may not 
Justify the effort for this particular task. 

(1) Assumption that s O. 

(2) Assumption that P^ in bubble is same as P.t> in emulsion. 

(3) No consideration of nucleation. 

(4) Neglect of flow contribution entering emulsion phase. 

This modeling is yet to be confirmed by experimental data. The 
experiments are underway. 
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APPENDIX 

Mathematical Model Program Listing 
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-For. IS SILIC0N_ „ 

REAL j 

C AAlRm. AN ARRAY, IS EQUIVALENT TO A I R( RPMIN+ 1 »0ELR-0. 5*DELR ) , A FUNCTION. 

C This is introduced to save the computer time by sto ring onc e co mputed numbers. 

C ALL' OTHER'S ARf 'o I MyLAR". 

COMMON AAIR(200 I , AABR < E-O ) , APO ( 200 1 • API 1 200 ) ,AP2 ( 200 ) 

COMMON AR^R(2C0),AAf;(20C) _ __ 

COMMON UO ,UMF , U f . OP »DT »H » CS » CS I L » CH2 , AL MF , aI F » Rh6s‘» RHOG . AMU 
COMMON EMF»ii.i.FO»Fl»F2»G,RE.OELP»RET»COP.UOR»NOR.TOH,AN»OB 

COMMON DE LR » UBR » UB . DE L T A » UF , ALPH A » US . . I » AKC_Ep S . UE tP SA .DSR . RC.Q 

■ COMMON AKOBA '.CABAV2 
COMMON jHMF 

COHMON_^|c:]^»5H_1.2_E_TA»a_p_* GAMMAE ,SHT»Alfp_BAy_tCA5AV 

COMMON EF » AKBCB » A<CEB » AKB EB » D. G AMMAC , GAMMAB , E T A , SH i AKDB » SC 
COMMON REP.ZAA.EBA.VA.AT.GRCWTH.PHIS .FR.AOR.Y.CAB 

COUPON AMW , RPMI N . RPMAXF , RPMAXB .TOT 

Dl .' ■'DyAMETER OF BED 

DP = PARTICLE DIAMETER 

EMF = MINIMUM FLUIDIZING VO I DAG E 

RHOG s GAS DENSITY 
RHOS»SOLID DENSITY 

AMU - VISCOSITY COEF F 

PHiS=SPHERICITY 

DElPMIN = PRESSURE LOSS ACROSS DISTRIBUTERS 

D “ DIFFUSI ON C OEFFICIENT (CM«*2/oEC) 

Y=L0GAKTTHMTC mean fraction ‘OF ■ fNERt ' gas' 
i.< = BED *EIGhT 

CAB=M£AN concentration^ OF SILANE IN THE BUBBLE (G R-M Q LE/CC) 

RPMIN = MlN'PARflCL'E SIZE 

RPMAXB = MAX PARTICLE SUE I\ LED 

RPMAXF = MAX P ARTICLE SUE IN FEED 

RRR‘ = PART'ICLE GROWTH RATE IN CM/' iC 
FI = OUTFLOW FROM BED 

PO = particle SIZE DISTRIBUTION - FEEDING 

PI = PARTICLE SIZE DISTRIBUTION - BED AND OVERFLOW 

P2 - PARTICLE SIZE DISTRIBUTION - ENTRAINMENT 

Fo = mas_s_flo'w rate - feeding 

F1‘ = MASS FLOW RATE - OVERFLOW 
F2 = MASS Flow RATE - ENTRAINMENT 

0 = GRAVITATIONAL CONSTANT 

arbitrarily assume The maximum particle size in bed. 

READ iOOC.UT.PHlS.RHCG.RHOS.AMJ.W.Fl . EMF, Y, CAB 

*, RPMAXF , RPMAX B. RPM IN, NOR >AMW,D,G ,5ELR _ 

tVERYIHING' IS EXP'RESSED IN CGS UNITS) NO CON'/FPSION FA'cTOR IS THUS NEEDED. 
PRINT ttda.DT.PHlS.RHOG.RH' S.AMU.W.Fl .E mE.Y.CAB 

RPMAXF. RPMAXB.RPMIN.NOR.AMW.O.G. DEEP __ 

sao FORMAT (Th T/7/- DT--,C1C.A.- phi £ = -, ElO.A,- RHOG=‘-.Fio. s‘, 

1- RHOS^-.ElO.i,,- AMJ^-.CIO.A,- 10. A,/- F1=-,E1C.A,- EMFa-.ElO 

Z.A.- Y--jE 10._A»- CAC=-.E:G.A,* RPVAXF = -,E1C.A»- RPMAXBr -,Fip.<n/ 

3 - RPM I N‘=- .‘ETO . A . -“NCR = - , E 1 C , A . A MW = - ,‘E 1 0 . A , - 'D r -T ,~E 1 0 , A - G « - 
A .E1^..A»- 0ELR=-,E10.A) 

DPM1N=RPMIN*2. 

II jso 

6 IIJ^IIJ+1 

IF(IIJ.EQ.I) GO TO 1-33 


PRINT COMPUTED QUANTITIES AFTER EACH CYCLE. 
RPMAXsRPMAXB 

NN=(RPMAX-RPMIN)/DELR .C.5 
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DO 151 1:^1»NN«5 

R JsRPMl N+ I *0CLR-0 . 5*OELR 

POjaAPOtl) 

l^lJ=APl(n 

P2J*AP2tI» 

RRRJ«ARRR( I ) 

AlRJ-AAIR(I) 

ABRJ^AABRU) 

PRINT 152tRJ»P2J»PlJ»R RRJtP0J»AIRj«ABRJ 

152 FORMAT <1H »-R“.P2.P1.0R/DT .Pb~J. AlPtABR=_,7E10.4» 
151 CONTINUE 

153 CONTINUE 

luOw format (8E10.4) 

IF(lIJ.NE.l) 00 TO 47 

C ARBITRARILY ASSUME THE MAXIMUM PARTICLE 
C ' SIZE IN BED MAYBE 1000 MICRONS * OPMAXS 
C 

C average particle SIZE)!) INITIAL _ 

ARP*0. 

NN*IRPMAXF-RPMIN)/DElR«^C.5 

DO 11 I=1»NN 

R =RPMIN+I*DElR-0.5*DELR 
ARP=ARP+P^' ( R) /r«delr 

11 CONTINUE 

RP=1./ARP 
PRINT 55» RP 

55 FORMAT! IHOt/- RP‘-.E10,4) 

GO TO 13 
C 

C AVERAGE PARTICLE SIZE) ) J ITERA’. .ON 

47 CONTINUE ' ' 

BRP*0« 

NNa<RPMAXB-RPMIN)/DELR »0« 5 

00 1.^ 1=1 »NN 

K =RPMIN+I*DElR-0.5«0ELR 
BRP=oRP^-P11R)/R*DElR 

12 CONTINUE 

RP=1./BRP 

PRI NT 55« RP 

c ■ 

13 CONTINUE 

C 

0P»RP»2. 

C AT * BED AREA 

ATaQT»»2»3«l4159/4 . 

C 

C COMPUTE UO. 

call TESTU 

C 

C DETERMINE THE TYPE OF FLUIDIZATION. 

CALL type 

C 

C DETERMINE THE DISTRIBUTOR PARAMETERS. 

CALL OISTRI 

C 

C COMPUTE TOM 

CAL L TDHS _ 

C 
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_C CHARACTERIZl BLiaOLES. 

CAUL 6UQBLE 
C 

^ CHARA CTERI ZE EMUL £-ION PHASE. 

CALL 'ErtULSl ' 

C 

_C CALCULATE GAS INTE R CHA NGE COEFFS. BETWEEN DUE9LE AND F^ULSlfi,N, 
CALL GASINT 
C 

C calculate HASS TRANSF ER COEFFS. 

CALL TRANSF 
C 

IF(llJ.NE.I) GO TO 121 _ _ _ _ _ 

NN«(RPMAXS-R?MIN)/OrLR+«.5 ' 

C 

OU 111 Ul.NN 
TF( JKL.EQ.21 uO to 211 
ARRR( I »=RRR(RPMIN-* I*DELR-C.5*0ELR) 

GO TO 111 

211 ARRRU )»RRR2(RPMIN+I*D£LR-0.5»OELR) 

111 CONTINUE 
DO 112 Ul.NN 

AAK (I)=AK (RPMIN+I*6ELR‘-0.5*DELR) 

112 CONTINUE 

DO 113 Ul.NN __ 

APU U'j^PO (RPMIN+!*DELR-0.5*Dt.LR) 

113 CONTINUE 

DO llA Ul.NN 

AAlR(n=AIR(RPMlN+l*6ELR-6V5*DELP)‘ ‘ 

114 continue 

DO 115 1 = 1»NN 

AABR ( I ) -ABR (RPM I N+ I »DELR-0 . 5*DELR 1 
113 continue 

121 CONTINUE 

C 

T0T=T0TAl(RPMAX3) 

r 

OiT 116 I = 1»NN 

API (I)=P1 (RPMIN+I*DELR-0.5*DELR 1 

116 CONTINUE 
C 

C BALANCE MASS FLOWS. 

CALL MASSFL _ 

C 

DO 117 I-ltNN 

AP2 m=P2 (RPMIN-* -I«DELR-0.5»DELR) 

117 CONTINUE 

C ONLY TWO ITERATION IS ENOUGH. 

IF ( I I J.LE.3) GO TO 6 

STOP 

C 

FUNCTION PC(Z3) 

c 'TNPT7T Particle oTzT' dT ser i t u t i on 

P0 = u. 

IF (Z3.GE.RPM1 N.AND.Z 3.LE.RP'^AXF ) PO - 1 . / ( RPF^AXF-RP M N ) 

return “ 

c 

c 


A-l» 
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C LUi aor^WJL.Tl.-;tL i :. ~ri C ALTLiLA lE IH£ 2Pfci^AT10atA;.. V E L0r . I7Y B A M 6 E « 

C ALMF - ritlOHl AT KF 

AcMr =ib/ J.i^i55»A./C:*»r/ ( 1 .-£MF)/ (RHOS-RHOr. ) /O 

utL.R=»< / S» in •? *«»» »i>T*»2 

C uMiL Vl K^*U DO.\-T K.\0» '£, SO Tmij is TRIAL AND ERROR -ETHOn. 

C i A^SOMt Rc.LT.iO 2 CCMPl-TE j;<F C TECIDl •JC(V'''’*UT I A CALCULATE RE 

C howeve r* Iw iniS k^RIICULAR CASE yO 2.,»JT 

RE * i 

{PASSsw 

L DELP = P RES SURE LC oS ACROSS 3E0 

L HL - REYNOLDS .'iC. 

V. u;-tF - Ml;, ri.UlDUING VELOCITY 

t UT = TEi^MlNAL V£ >.DC1TY ... 

11 IPASS * IPASS^l 

OMF = yP»*2«tRhOS-RMOu)*C/165C./AMU 

uT = Ji«(_RHpS-RHOC)»DPMIN**2/12./AMU 

Ir ) JMF = saRT(Dr*jFHOS-RHOG:*n, 24.5/RHOGj 

IFIRE. CE.2S.J PRINT 51 

C USE. ANY U C. JUST PLACE THE.?15tR£0 ONE IN THE LAST POSITION, 

ULs6,5»UT 

U0-.25*UT 

UDf.l_.A*UI._. - - 

uO=2.AUf 

UMF2=2.*UMF 

UO-2.5*_yT .. 

UO*19.'23 ' 

1F(Uw»lT.JHF2) oCsUKF2 

RL=pP»RhpG*oL/AMU 

1 * I IPASo.EO. 1 ) so TO 11 
51 format RE IS lARSc>\ T;^A,'4 20—) 

IFIoMF.LE.uT) GC ’0 32 

21 PRINT 31 

31 FORMAT IlHo,- UT IS Sr^ALLER Tt'.AN UMF-) 

32 CONTINUE 

PRIN 61.0MF »UT tREtAuMF »DELP,J0 

61 hORHATIlMS,- LHF ,uT »R£ »AL.Vr ,DELP »U0 = - »6E 1 3. A ) 

A1 RETURN 

C 

C 

subKOuT ine Type _ 

L FR - FROUDE NO 
L REPKr » RL AT MF CONDI Ti ON 

C mNUMuER CRlT^'^iC.'. OF TYPE OF FLUlO ! :/• TION 

C TnlS SUbROuTINE'iS TO CHEC< THE DECREE CT FLUIDIZATION IN PED 
FR = UMr**2/DP/G 

REPMF a oP»UMF»RKCG/AMu 

ANuMcd-i f\*' REPMF •(Ri;OS‘*RhCG) /.\H0G p AL''? »’D T 
i F I ANUMdE .SE, iw w • < PRiNTl 

IF I ANyHbE . L T • i 00 • ) GO TO 5 

1 FORMAT attO»////- 0U3bLUNG FLJI DI2AT ICN-) 

UTEST a 2*UMF 

IF ( Uu.GE.uTEST ) PRINT 2 . .. 

2 FORMA! THIS IS VISCROUS-) 

I F i Uw. L T. OILS r ) PRINT 3 

;► FukKAT (1H0»- THl- IE NOT VIGOROUS-) 

5 COMINwc. 

I F I ANUi’iOL .L T • 10 0 • i PKlNl H 
A KjhHAT PAR'iCULAiC v LoIO I ZAT i ON— ) 

1 F ( ANyi'lO. CE • 1 • ) v!^u“i 
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• 1 h I . .. 

kHlukU 

C 

c _ .. 

SUbROLT ii\£ 0I3 Ti7I 

THIS Sul>RC.u; INw IS rc CALCULATE Tht C I S T*< I l^UTOR 

velrat ^ VLLCCiiv r.Ario 

AOR = CRIFICE AREA 
CDP = GRIT ICE COEPriCIEM 

NOK - HuHbLR OF OR I ICE ... 

UOK ^ OKI FILE velocity 

DElPi ^ delp»o.i 

DclP^=3 . 

utLF^ io LOLIvAlENT to 35 CA H20 
REFER 10 P67 Or RUXII» EQ 3A. THE THIRD TERM 

1 5 imcOELELTED ^^ERE 

OELPD^uElPI 

IF tuELP2«oE.0LLPi: 0LLrD=DELP2 

iFtDuLP2*GE.0ELPl.> ^RJNT 11 . . . 

ii FuRMAi (Y hw#///- *E CAn INCREASE HEIGHT OF TH> SEO wlfHOHT PENALIZING THE 
♦DESIGN ulSTKIBJTCR PRESSURE LOSS*) 

KtT - DH*KH0G*0O/A.Xu 

REFER TO FIG 12 OF" Pda 

iHET*GE*30«) CD?~S*6 

1 F I Rw: . LI . 3v . ) COP = C . A* ( AL0G( RET ) -AL9G ( 2* ) .) / ( ALOC-f 30# ) -ALC G( 2> ) ) 

UOK = COP* SORT i 2# -DElPD/RHOG) 

VElRAT-Uu/oOR _ 

AOk - VELrvAT/NCK 
PRINT 2 A ♦ RL T ♦ COP »UwI\ » \/EL RAl »AOrj 

2x format i •j-t3£: fCDP, jO.^ , V£LRa:» = - * 5 E IO#A 

RETURN 


subroutine TdhS 

In-S cUbROUTINE IS TC CALCULATE T?H 
REFER lo P9A 

lOn - IKANSPCRT DISENGAGING HEIGHT 

uwr)K'i Sv^bRwUTlNE 

RETURN 


SUbivOw 1 w>*.^bbLE I 

UF oFWA;0 V'lLCCITV cf gas I\ emulsion AT Mr 

uElTA = VOLUME FRACTICN :F SUGGSES 
KB = RADIUS OF BUBBLE 

TrtiS SuBROs^TINl COMPUTES EUESlE PROPERTIES# 

RErCiN 10 LriAP A 

c> b * * Ai*l^ f L (V i‘ L wL oL fL 

an Soi^uLL i RL«« jL\Cy 

UbR * subtle RI Sl,‘,S V:.^CC I TY 

Ub = ABSOLUTE RISE .ELCCITY 2. BELTS 

H = height BETWEEN TWO SUCCESSIVE PUPBl^C 

sc-gcmmidt no. 

rUr% simplicity OF THE CALCULATlCNf BASED ON TMF ! T | 

Fiu 2u*Pl29l LCT -S assume THAT THE BUBPLF '“PEQUENCY IS APPPOX T MA T'ei Y 
4 regardless of the size of the bed ^OR the most part cf the BED# 

AN - H. 

DB* 1 • 5 / AN^ ( vw - uxr ) 
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- . . . 

1 FORMAT (IhC.//- CJ8BLE ^UE. FRE3UENCY --•7E10.4) 

OQK - 2^«26 * DB»» (w.Sl 

C the bubble size cannot be larger than max stable SIZE* EO.o.PlZl 
IF( UBR.Gt.UT) UBTRsUT 
OB«UBR**2/22,26*»2 

RB»DB/2. 

UB»UO-UMF+UBR 
H > UB/AN 

UF « UMF/EMF 

DELTA=2./3*»DB/H 

PRINT 2»UBR»UB*H»UF.0ELTAtDB 

2 FORMATUHG*/- UBR>UB*H . uF* DELTA. B=-»6£10.A) 

RETURN 

C 

C 

SU3R00TINE EMULSI 
C EMULSION phase CHARACTERISTICS 

C REFER TO CHAP 4 ANO 5 _ 

C RC s RADIUS Or CLOUD 

c Q * Flow rate of gas into and out of a bubdle 

C (P115) 

C ALF - HEIGHT Or BJEolING BED CP131) 

C EF * VOID FRACTION IN SUSOLING SE01P131) 

C UE = UPWARD velocity OF GAS IN EMULSION 

C ALPHA = RATIO OF WAAl' VOLUME TO BUBBLE VOLUME 
C US = MEAN DOWNWARD VELOCITY OF SOuIDS IN EMULSION 

C AKC£dS(P 157) - COEFF OF SOLID INTERChANGt DZTWLEN 
C ClOUO-WAaE region and EMULSION PHASE 
C DSA = AXIAL OISPLRSION COEFF OF SCLIOS 

C DSR = radial disp ersion COEFF Or SOL ICS 

RC » 0*5»DB«( (U8R+2'.»UF )7(U3R-jF ) )*“^S.3333 
C OBMAa would not be defined IP122 -"IG.14) 

a = 3.»UF*EMF » 3.i O*Dti»»2 - 0,23 

ALF = Ali»ii / (i,“DELTA) 

EF = l.-i 1, -DELIA ))Ul,-E.MFi 
PRINT 1» RC» a* ALF * EF 

1 FORMAT (IHO./- RC,G,ALF,EF.=-»^E10.A» 

C FOR SIMPlICITY»BASED ON FIS,9,P150, LET-S ASSUME THAT ALPHA IS CONSTANT 
ALPHA=U»A _ _ 

US = ALPHA*DELfA»u-/li.-CELTA-DELTA»AEPHA) 

OE=UF- J j 

J = ALPHA>DEwTA*>/d«RH OSM 1»-EMF ) 

AACEBS - 3,*« 1,-£MF)»UMF7 C l.-DELTA)/EMF/Da 
0 SA=AlPHA»«2»EMF*0E*(U0-UMF)»«2/3./DELTA/UMF 

_ DSR = 3»/1o. *DELTA/(1.-DELTA)»'JKF*DS/LMF 

PRINT 2* ALPHA. US, J,AACEBS»UL»DSA,CSR 

2 FORMAT ( 1h^-./— ALPHA. Uw.j. A A CEoS .Uw.DSA.DSR — ".*^u1U,A) 

_ return 

c '' 

c 

subrouti ne g asint 

C GAS INTERCHANGE OETwEE.'. bJBBLE AND CMULISON, 

C WILL BE applied TO ThL CASE CF VIGOROUS 
0 BUB BLING BEDS EITHE R U0.GE.2UMF OR 
'C 'ud.GE.5UF 

C AABCB = interchange COEFFICIENT BETWEEN 
C BUBBLE AND CLOi^D 

C AACEfa - dETwEEN CLOoO A,N0 EMULSION 
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Cl a a£TM£EK BUBBLE AND rttULitrOM 

C THCSe COCFFS* ARE A&SVMEO TO BE INOERENCENT OF IKOXVlDUAU PAST10.E SIZE 
AIC8Cd*A*S»UMF/Od»5*8S» C 0»»2«G/0Q« »5 » • « { C»25 > 

AfcCEB ■ »«7> * (EH F»DOUB/SB*»3»»»ta.St 
AK6C&*l«/(l./AKBCBn«/AiCCEVr ' 

PRINT If AKSCBfAKC.LtfA<B£a 
1 format UH w,/- AKBCSfA KCEQfAlca£B--.2£ 10.41 

return 




C 

C '■ 

&U8RCUT1NE TRANSF 

C OBTAIN DIFFUSION COEFFICIENT ON TO S1«.1C0N SOLIDS 
ETA»1. 

C Cf^fECTIVEfiEsS FACTOR 3 

C ETA IS TAKEN TQ BE !• SINCE THE PROCESS IS CIFFUSION C(HITR0LL£0 
C AND THUS THE SILA NE CONCENTRATION CM SOLID PARTlCLESfCAS* IS NO! 
X IKPORtANT IN OUR CASE* 

C SAN4Ad> VOLUNE FRACTION OF S(».1DS IN 6u3BLE 

^ ■6ANHAB iS AR&IlRARILY CHOSENIPaOJi - . 

GAI1NA0 n_no3 

C CANNACa \»OLUNE FRACTION OF SOLIDS IN CL(AN> 

6ANNAC« C l.-E.’y i »C 3.*UNf /ENF/ IUQR» tdtf^-FKF i ^ALPHA i 

C GAMKA£« VOLIME FRACTION OF SOLIDS IN EMULSION ^ 

GANMAE*(1*-EHF) •{ 1 ••DELTA) /DELTA' (GAMMAB-»^CAMT'.AC) 

C SC _* SCHHIDT 

SC*AMU/RH06/D ‘ 

C RETT « REYNOLDS NO* NITH UT 

RETT«DP*UT* R HOG/AN U 

Z REMF * REYNOLDS NO. WITH UHF 

REMFs DP«UKF»RHOG/ANU 

C EETA a A VARIAOLE INTRODUCED IN NEL SON S GALLJiiiAr. kilTH vfl 

EETA«(l»/(l»'EF)«»(l«/3*)'l.)»D*S»RE*«tO*5)«SC*»ll*/S*) 

C ZEMF = SAME VARIABLE WITH UMF 

2ENF=Il»/Il»'EF)»«Cl. /3.l'l.l»Q.3«R£MF«»(Q.Si«SC»»Ily/^.) 

C SHMF = SHERWOOD NO. -ilTH NiF CONDITIONS. TO BE USED IN i-MULSlOW ' 
C PHASE In BUBBLING BED 


C 


C 

£ 

c 

c 

T 

c 


C 


SHMF»t2.»ZEMF^-t2.»EEMF**2*U. -EF)*« ( 1 ./3. )/t l.- l I.-EF M »a./5.l l 

♦••2-2.I*TAN)IIZEHF|)/lZEMF/<l.-< l.-£Fj»*Il./3.n-TANHCZEMF) J 
SHI « SHERWOOD NO. WITH UO. TO BE USED IN PARTICULATE uED. 

SHI ■(2.»ZETA*<2.*2£TA»«2 »(1.-EF> 1./3. ) / 1 1.-< l. -EF) «»( T./l.n 
♦•*2-2. >»tANHliETA) )/|ZETA/tl.-( I.Vef |»»c i.y3. ) )-TA.’iriCZElA » ) 


SHT * SHERWOOD NO. WITH L'T, TO BE USED 
SHI»2.0-»0.6»SC»»(l./3.)»SQ?niRt:TT} 


FOR SOLIDS IN bubbles. 


AK0BA2« AVERAGE MASS TRANSFER COEFF. FOR PARTICULATE BED. 
AK06A2«SHI*0/0P 


CABAV2»AyERAGE CONCENTRA TION OF SILANE IN PARTIC ULATE BED . 

CABAV2«U0»CAB* C I .-EXP ( -AKDBA2*ALF/U0 ) ) /AKDBA2/ALF 

COMPUTE AVERAGE SILANL CONCENTRATION IN BED. REFER TO KY WRITEUP. 

DP IS THE AVERAGE PARTICLE SIZE. 

80 * AN INTERMEDIATE VARIABLE. 

B0«6 • »0»SHMF/PHI S/0P«»2 

AKDBAV* AVERAGE MASS TRAN SFER COEFF. FOR BUBBLING «££* 

AKb8AV«GAMMABiaD«SHT/SHMF+r./7l.7AK3CB*l./iGAHNAC*B0+li 
*/t l./AKCEBf>l./GAMMAE/BD) ) ) 

CABAV »AVi:RAGE CONCENTRAT ION OF' SILANE IN uUDBLING JED. 
CA8AV*Ua*CAB»( 1 .-tAPl 'AKOBAVvALF'/UB i V / AkOBAV/ALF 
PRINT S»CABAVfCABAV2 
3 format ( iHOt- CABAV fCABAv2—f2E10.4) 

“RrroRU 


rr^ nnnnnn r>rv ; r%n nn 
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FUKCnOH AIRIA) 

AIR ^ INTEGRAL CQ.32. P339 
10 COMPUTE THE TMTEGRAt KR.Rtl 
»M«tA-RPMli4l/0Et.R^w.:» 

A1KLOO*<<« 

OO 1 

R9sRPKI a* I •OELK-0. 5«KCR 
AlRLO&>AlRLOG«tFl/lk-frAAKI 1 } l/ARRR 1 1 I«0ELR 

i COKTIWUE _ 

AIR-EXPT-AIRLOG) 

RETURN 


FIMiCTlON 4^1211 
COMPU TE AN INTERGRAC 

ABR - INTEGRAL ShOUN IN EQ.36. P340 
ABR^C* 


00 51 I>1*NN 

R1 »RmiN-i-l •DELR>0.5«0 ClR 

ABR« AaR*AP0m/Rl»«3 /AAlRl UfDELR 


51 continue 
RETURN 



' FUNCTION TO'TAlCZAI 
COM>UTE AN INTEGRAL W/FO IEQ.37.P340I 
TOTAL-0* 

NN» UA-RPMINI /OELR+0.5 

DO 1 I >1*NN 

R6« RPMIW » 1 c OELR-0 • 5«0ELR 

T0TAL«TCTAL4-R8«*3/ARRR( I I«AAIR: I )«0ELR»AABR( I ) 

1 continue 

RETURN 


SUdROUT InE MAS3FL 

■ OElR =^PARriCL£ 5“I2E INCRE>-I£NT USED FOR 
COMPUTING the integrals 
FO< hi/TCT 

COMPUTE Fw*f2*P1»P2 

GROWTH « particle GROWTH RA»E IN 3£0<C3.23.r3i7l 

GHOwTH ap _ 

i RPMAXO-i^M I N )70 ClRVc • S 
OO 1 1^1*NN 

R6«RPKI t\* I aOELR'S • 5*0 ElR 
0R0wTtUGR0WTH^0CLn»3.*W»Ari( 1)*ARRR( 1 )/R6 

1 CONTINUE 

^ GRO*Th - Fl+fC 
PRINT 2* F0iFl*GR0WTH*F2 

2 FORMAT (IHL.//. F0»Fl,GR0i^TH»F2>-*E£I.d.A> 
RETURN 


C 


function .3RR(X1» 

TO compute 6R/0T for BU36LIN6 BED 

X»2*»Xl 

RRRa(DeLiA*X«PHlS/e./(l,-cri}MCABAV«AM*/RHCS)«(GAMKAe«S«*D»ShT/ 
♦PHlS/X»»2+l»/(l*/AX£CB+l./lGA'?MAC<6.*0*SiiMF /PHIS/X«»2 + 
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OMQJNAD PAOBJB 

OR POOR 


. Vil,/AKCL5'l«/(0A.'tMA£:»6. 0»SH'‘K/Fn:o/X*’»25 J : J ) 

k&IURN 

c 

fUNCilO<« RRR^iAl 

C 10 CoMFoIE Oiv/CT PAnliCoLAfE uLO 

RtUlRH 

C 

FJNCriOA AX(A) 

C Inli COMPuic- t'HL ELUTRIAIION COErrICIEMT 

V, FiiOH IRANiFE RRA=^AA/iP tiSA^DF. 

V Xii - REFER 10 F314>» flG 13 

C _ Ai^ - LLOTRIATIOM CONSTAi«I AS A rONCTiCN OF RADIOS 
C ' ASI - tUC-UT)/SCKUGl»(tJI»RHOC/A.MUI*»C«:25 
C •( IRhOS-ahO&i/RHCGJ’H'i.lS * 0P»*C,225 

C = VA*OF**0,225 

C hO E4.oTRiATlC.'4 IF UO.UE.oT 

IFtOw.LE.un GO T(. 1 

VA =10o-uX»^^Rr‘G>»(OT*KHCG/AM0; 3 «u,7£S 
»•! (RflOS-RHCOl/RHOGi^ilalS 

AA «EXP*ALOGIO«15>*<AECG£17.» -AcCG(C.lS; 1«{ ALCC£ VA* * 2« sA« »» 
AO • 225 } — ALoGC 100CC«# #/( ALCGf 10 CCOOO • « — ALOC ( 10000* i « i*Af«00 

• viT J *AT /ini 

GO TO 2 

1 COTiTINwL 
AA=« 

2 CONTINUE 
KETURT<i 

C 


FUi%CTi\^i« rlCZoi 
C cG«3S*P3AO 

Pls26*»3 / ARaRI ; ) SAAIRC I : K AASR t 1 | /TCT 

return 


Function pzizioi 
C EO«7»'P329 

P2^AAA( 1 |•M/F2•AP1 ( 1 ) 

RETURN 

END 


-HAP 
-XU I 









£♦01 

• 66 

£+00 •1236E-03 £+0i 

. 3ASAE -03 

.196 E4-06 

•9 L+03 

• 3 £-»vu 

• A 

£♦00 

• 3 

£“03 •TS £“02 *9 £-01 

•A5 E-02 

.1 C-*-03 

• 2iiC0L*C,Z 

• 3607C + 0a. 

• 96 
-Fils 

£♦03 

• 5 

£“03 
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